Parallel Adaptive Finite Element Algorithms for Solving the Coupled Electro-diffusion Equations

نویسندگان

  • Yan Xie
  • Jie Cheng
  • Linbo Zhang
چکیده

In this paper we present parallel adaptive finite element algorithms for solving the 3D electro-diffusion equations such as the Poisson-Nernst-Planck equations and the size-modified Poisson-Nernst-Planck equations in simulations of biomolecular systems in ionic liquid. A set of transformation methods based on the generalized Slotboom variables is used to solve the coupled equations. Calculations of the diffusion-reaction rate coefficients, electrostatic potential and ion concentrations for various systems verify the method’s validity and stability. The iterations between the Poisson equation and the NernstPlanck equations in the primitive method and in the transformation method are compared to illustrate how the new method accelerates the convergence of the solution. To speed up the convergence, we introduce the DIIS (direct inversion of the iterative subspace) method including Simple Mixing and Anderson Mixing as under-relaxation techniques, the effectiveness of which on acceleration is shown by numerical tests. It is worth noting that the primitive method fails to solve the size-modified Poisson-Nernst-Planck equations for real protein systems but the transformation method succeeds in the simulations of the ACh-AChE reaction system and the DNA fragment. To improve the accuracy of the solution, we introduce high order elements and mesh adaptation based on an a posteriori error estimator. Numerical results indicate that our mesh adaptation process leads to quasi-optimal convergence. We implement our algorithms using the parallel adaptive finite element package PHG [53] and high parallel efficiency is obtained.

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Fast Finite Element Method Using Multi-Step Mesh Process

This paper introduces a new method for accelerating current sluggish FEM and improving memory demand in FEM problems with high node resolution or bulky structures. Like most of the numerical methods, FEM results to a matrix equation which normally has huge dimension. Breaking the main matrix equation into several smaller size matrices, the solving procedure can be accelerated. For implementing ...

متن کامل

Parallel Finite Element Method for Coupled Chloride Moisture Diffusion in Concrete

Penetration of chloride ions into concrete and diffusion of moisture in concrete are important factors responsible for the corrosion of steel in concrete. The two diffusion processes are coupled. This paper deals with the analysis and simulation of coupled chloride penetration and moisture diffusion in concrete. Of particular interest is the parallel programming in finite element method for sol...

متن کامل

Theoretical Formulations for Finite Element Models of Functionally Graded Beams with Piezoelectric Layers

In this paper an overview of functionally graded materials and constitutive relations of electro elasticity for three-dimensional deformable  solids is presented, and  governing equations of the Bernoulli–Euler and Timoshenko beam theories which account for through-thickness power-law variation of a two-constituent material and piezoelectric layers are developed  using the  principle  of virtua...

متن کامل

A Closed-Form Solution for Two-Dimensional Diffusion Equation Using Crank-Nicolson Finite Difference Method

In this paper a finite difference method for solving 2-dimensional diffusion equation is presented. The method employs Crank-Nicolson scheme to improve finite difference formulation and its convergence and stability. The obtained solution will be a recursive formula in each step of which a system of linear equations should be solved. Given the specific form of obtained matrices, rather than sol...

متن کامل

Parallel Adaptive Simulation of Coupled Incompressible Viscous Flow and Advective-Diffusive Transport Using Stabilized FEM Formulation

In this work we study coupled incompressible viscous flow and advective-diffusive transport of a scalar. Both the Navier-Stokes and transport equations are solved using an Eulerian approach. The SUPG/PSPG stabilized finite element formulation is applied for the governing equations. The implementation is held using the libMEsh finite element library which provides support for parallel adaptive m...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

عنوان ژورنال:

دوره   شماره 

صفحات  -

تاریخ انتشار 2013